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Introduction 


This is the final report of the NASA Lewis SBIR Phase II Contract Num- 
ber NAS3-25785, Multigrid Solution of Internal Flows Using Unstructured 
Solution Adaptive Meshes. The objective of this project, as described in the 
Statement of Work, is as follows 

The objective of this SBIR Phase II project is to develop and 
deliver to NASA a general three-dimensional Navier-Stokes code 
using unstructured solution-adaptive meshes for accuracy and 
multigrid techniques for convergence acceleration. The code will 
primarily be applied, but not necessarily limited, to high speed 
internal flows in turbomachinery. 

The objective of the SBIR program, as described in the NASA Lewis Re- 
search Center Annual Report at the time this project began, is as follows 

The Small Business Innovations Research (SBIR) Program is 
a congressionally mandated program aimed at the commercial- 
ization of innovative concepts by small business. These innova- 
tions are in general areas of research based on the needs of the 
U.S. Government. 

This contract has successfully achieved the objectives of both this contract 
and of the SBIR program. On June 1, 1992, Fluent Inc. introduced a new 
Computational Fluid Dynamics package, called Rampant, into the industrial 
and research markets. Rampant is an unstructured, solution-adaptive code 
that can solve inviscid, laminar, and turbulent; compressible and incompress- 
ible flows. The initial development stages of Rampant were made possible by 
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Chapter 1 

Overview 


The licensed copy of the Rampant package provided with the contract deliv- 
erable includes a copy of the Rampant User’s Guide. The User’s Guide 
provides a comprehensive description of Rampant’s capabilities and opera- 
tion. This report does not, therefore, include a detailed description of the 
mechanics of operating the program, but instead concentrates on describing 
the numerical techniques and test results that are pertinent to the specifica- 
tions of this contract. 


1.1 Program Capabilities 

Rampant has the following general capabilities: 

Grid Generation 

• 2D unstructured triangular mesh generation. 

• 3D unstructured triangular surface mesh generation. 

• 3D unstructured tetrahedral volume mesh generation. 

• Can read triangular surface meshes from a variety of CAD packages, 
and can then mesh the interior with tetrahedra. 

• Can read tetrahedral volume meshes from a variety of CAD packages. 
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• Arbitrary slices. Planar, spherical, cylindrical, etc. 

• Iso-surface extraction. Extract surfaces of constant value. 

• Iso-line extraction. Extract lines of constant value on a given surface. 

• Surface integration. Pressure and viscous forces. Normal and shear 
forces. Lift and drag coefficients. Heat transfer. 

• XY plotter. Residual histories, surface plots, over-plot data from file. 

• Generate Tecplot and Cray MPGS format visualization files. 

1.2 The Solution Process 

The process of solving a flow problem with Rampant consists of the following 
primary steps 

1. Generate a geometric surface description 

2. Generate a surface mesh. 

3. Generate a volume mesh. 

4. Setup the problem to be solved. This includes setting boundary con- 
ditions, flow parameters, and solution parameters. 

5. Solve the flow. 

6. Analyze the flow field. 

7. Adapt the grid. 

These processes are described in the following chapters. 
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• Turbulent Kinetic Energy 


i JJJ pk dV + JJ pk VidAi - # S S ^ = JJJ (P ~ pe) dV(2A) 

• Turbulent Dissipation 

I III p(dv +jj^ iiA ' - JJ ~M = 

JJJ -^(ClP - C 2 pc)dv ( 2 . 5 ) 


Euler Equations 

For inviscid flow, the momentum and energy equations reduce to 


• Mass 

lJJJ p Vjdv+ JJ pvj VidAi + JJ pSij dAi = JJJ pgj dV (2.6) 

• Momentum 

JJJ pEdV + JJ pE VidAi + JJ pvi dAi = 0 (2.7) 

2.1.2 Rotating Cartesian Formulation 

For rotating domains, the equations are formulated in terms of absolute 
velocity in a frame of reference rotating with angular velocity Q about the 
origin. The formulation is valid for Q ~ S2(t). The relative velocity, w, is 
defined as 

w = v — ft X x (2-8) 

Euler Equations 

• Mass 

diJJJ pdV+ JJ pWidAi = 0 (29) 
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where 


$ = Vv : (Vt> + (Vt>) T ) - 2/3 (V-t>) 2 

= VvjiCVvji + Vt»y) - 2/3(Vi?,,) 2 (2.15) 

= Vr*i(V«jci + VtJjar) + 

4- Vi?i y ) + 

V^(Vr 2i + Vi; tz ) - 2/3 (Vt,,) 2 
and for a Cartesian coordinate system, Vu is given simply by 

Vn t j = dvi/dxj (2.16) 


The k-e model parameters are 


fie = 

/^laminar 4" /^f 

Pt = 

pC^/e 

°k - 

1.0 

cr c - 

1.3 

Ci = 

1.44 

c 2 = 

1.92 

For laminar flow, 


pe 

= /^laminar 


fit = 0 


2.1.4 Equation of State 

For compressible flow, the ideal gas equation of state is used to relate pressure 
and temperature to the conserved variables. 

p = pRT (2.17) 


where 


R • — Cp Cy 
7 = Cp/c v 
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becomes 



pv • dS — 


X JJ p v • 6(4 

faces 


X // /’»• dA 

faces 


Z 

faces 


( 2 . 21 ) 


where the face-averaged value of velocity is given by 

pv 


v = 


and the components of the area vector, A , 


-ii 


dA 


( 2 . 22 ) 


( 2 . 23 ) 


are the area projections of the face. 


The remaining integrals are discretized analogously, and the Euler equations 
reduce to the following system of ordinary differential equations 


where 


j t (VW)+ £ F A = 0 
faces 


' p 


' p V 

fox 


pv x v + pi 

P v y 

, F = 

pv y v + pj 

pv z 


pv z v + pk 

. P e . 


pe v + pv 


( 2 . 24 ) 


2.2.2 Added Dissipation 


The finite volume scheme described above is analogous to to a central dif- 
ference scheme on a structured grid. The use of central differences ensures 
that the scheme is second-order accurate for “smooth” meshes, but allows 
the solution to decouple at odd-even points. 
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where K 2 and K 4 are constants, for example, 2 and 0.05. 

The application of the first-order dissipation is governed by the TVD-like 
switch of Turkel and Swanson [TSVW91] which is sensitive to second differ- 
ences of pressure 


2^(p+-p~) 

~ faces /n 01 

|(p + -p_)|+ J](p + +p_) 

faces faces 

This switch also turns off the fourth difference dissipation in the neighbor- 
hood of a shock where it can be destabilizing. 

For multigrid calculations, the blended second /fourth- difference dissipation 
is replaced on the coarse grid levels with a simpler second-difference dis- 
sipation. The dissipation is only first-order accurate, but accuracy is not 
necessary on the coarse grids and the second- difference dissipation is more 
stable and cheaper to compute. It is equivalent to setting 


£2 — ft 2 c (2.32) 

£4 = 0 (2.33) 


where H 2 c is typically 0.5. 

The face-based dissipation operator, D, results in added flux-like terms that 
are treated in the same way as the physical Euler fluxes. This results in a 
discretized equation of the form 

^-(VW) + £ (F • A - D W) = 0 (2.34) 

CLL r 

faces 

Since the facial dissipation operator, D, is a flux-like term, it is set to zero 
along external boundary faces to maintain global conservation. 

2.2.3 Divergence Theorem 

Extensive use is made of a form of the divergence theorem (or Green’s 
Lemma) in calculating derivatives and metric quantities. The divergence 
theorem is given by 

JJJ'V-FdV = (2.35) 
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A standard set of coefficients for a four-stage scheme is 1/4, 1/3, 1/2, and 1. 

This type of multistage scheme is, at best, only second-order accurate in 
time. For steady solutions, however, it is usually preferred over more accurate 
Runge-Kutta schemes because it requires the solution to be stored at only 
two time levels. Higher-order temporal accuracy is not required since the 
scheme is used simply as a means of advancing the solution towards its final 
value. A classical higher-order Runge-Kutta scheme could be implemented 
if the increase in time accuracy is deemed to warrant the increase in storage 
requirement. The classical fourth-order Runge-Kutta method is given by 


w (0) = w t 

IV (1) = W (0) - l/2AtRlV (0) 

= w (0) - l/2AiRVV (1) 

W (3) = U (0) - A(RVV (2) 

W W = W <0) - l/6At(RW^ 0) + 2IUV (,) + 2RW' (2) + Rff 1 ’ 1 ) 
Wt+At = 


2.2.5 Implicit Residual Smoothing 

The maximum time step can be further increased by increasing the support of 
the scheme through implicit averaging of the residuals with their neighbors. 
The residuals are filtered through a Laplacian smoothing operator. 


+ ( 2 ' 38 > 


The resulting set of equations can be solved by using Jacobi iteration 




Ri + eER'i l ' m ~ 1) 

i + fEi 


(2.39) 


Two Jacobi iterations are usually sufficient to allow doubling the time step 
with a value of e — 0.5. 



interpolation and relaxation operators. A sample coarse-grid hierarchy is 
shown in Figure 3. 

The problem that remains is the solution of the governing equations on the 
coarse grids composed of irregularly-shaped cells. Since the discretized equa- 
tions presented in the previous sections place no restrictions on the number 
of faces a cell is composed of, this actually poses little additional problem. 
There is a loss of accuracy when the finite volume scheme is used on the 
irregular coarse-grid cells. However, the accuracy of the multigrid solution 
is determined solely by the finest grid and is therefore not affected by the 
coarse-grid discretization. Also, a simple, more stable second-difference dis- 
sipation can be used on the coarse grids since they do not affect the accuracy 
of the fine grid solution. Higher-order terms such as viscous stresses and tur- 
bulence coefficients are “frozen” at their fine grid values and hence need not 
be calculated on the coarse grids. 

One further note of interest is that although the coarse-grid cells look very 
irregular, the discretization cannot “see” the jaggedness in the cell faces. The 
discretization uses only the area projections of the cell faces and therefore 
each group of “jagged” cell faces separating two cells is equivalent to a single 
straight line connecting the endpoints of the jagged segment. Although not 
currently implemented, this fact can be used to decrease the cost of the 
coarse-grid iterations by reducing each group of face flux calculations to a 
single flux calculation using the group’s net area projection. 

2.3.2 Multigrid Cycle 

A multigrid cycle can be defined as a recursive procedure that is applied at 
each grid level as it moves through the grid hierarchy. The traversal of the 
hierarchy is governed by three parameters, /? 2 , and /? 3 , as follows: 

1. fii smoothings , or relaxation sweeps, are performed at the current grid 
level to reduce the high wave number components of the error. In the 
present case, the smoothing procedure is the multistage scheme with 
local time stepping, residual smoothing, etc. The high wave number 
components of error should be reduced to the point that the remain- 
ing error is expressible on the next coarser mesh without significant 
aliasing. 
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In the current implementation, interpolation of corrections is per- 
formed simply by setting the correction in a fine-grid cell equal to 
the correction in its associated coarse-grid cell. Following this, the 
interpolated corrections are smoothed with one pass of a Laplacian 
smoother (as for implicit residual smoothing) to reduce its high wave 
number components. After being smoothed, the corrections are then 
added to the fine-grid solution. 

This flow of control can be expressed in psuedo-code as 

multigrid_cycle(n) 

{ 

for i=l to fii 

multistage Jteration(n) ; 

if (n ^ coarsest -grid Jevel) 

{ 

compute_residuals(n); 
restrict_solution(n,n-t- 1 ) ; 
restrict_residuals(n,n+l); 
for i— 1 to /?2 

multigrid_cycle(n+l); 
interpolate_correction(n+l ,n) ; 
for i=l to /?3 

multistage _iteration(n) ; 

} 

} 

The parameter that each of the functions calls is the grid level for which the 
result is computed, beginning with 0 on the finest grid and increasing by 1 
for each successively coarser grid. The highest level of control consists simply 
of performing multigrid cycles at the finest grid level until some convergence 
criterion is met. 


2.4 Boundary Conditions 

A variety of boundary conditions are available in Rampant. Boundary con- 
ditions are specified for groups of boundary faces that axe identified during 
the grid generation phases. 
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For rotating domains, pressure is specified at the “hub”, and a radial pressure 
distribution is computed using the radial equilibrium condition 

fe = eYl (2.45) 

or r 

where V$ is averaged circumferentially. Circumferential averaging is per- 
formed by dividing the interval from the minimum radius (hub) to the max- 
imum radius (shroud) into some number (e.g., 64) of evenly spaced “bins”. 
Each boundary face is then associated with the bin that corresponds to its 
radial position, and averaged values are computed for each bin. A pressure 
for each bin can then be computed by marching out from the hub and ap- 
plying the radial equilibrium condition to each bin. The pressure applied at 
each boundary face is then simply taken as the pressure in the associated 
bin. 


Supersonic Inlet 

All quantities are specified. 


Supersonic Outlet 

All quantities are extrapolated. 

2.4.2 Far-Field Boundary Conditions 

For external compressible flow such as that about aircraft configurations, 
far-field (or free-stream) boundary conditions are usually used. Free stream 
Mach number, static pressure, static density, and flow direction are pre- 
scribed and boundary values are computed in the standard way using Riem- 
man invariants (see, for example, [Jam83]). 
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• maximization of the minimum angles in the grid cells, the produc- 
tion of the most equilateral mesh for a given distribution of vertices 
(equiangular property); 

• simple criteria used to drive the triangulation process, either the max- 
min angle, the circle (2D) or sphere (3D), or the Voronoi neighbor 
criterion; 

• advantageous environment for mesh refinement provided by specific 
implementations of the technique; 

• flexibility gained by decoupling the grid point generation from the 
process of establishing connections between the points. 

The present approach for generating the interior mesh composed of either 
triangular (2D) or tetrahedral (3D) cells was proposed independently by 
Bowyer [Bow81] and Watson [Wat81]. The scheme performs a nucleus of 
operations, casting vertices into an existing Delaunay triangulation and us- 
ing efficient search procedures and the circle (2D) or sphere (3D) criteria to 
reconstruct the grid. The additional benefits of this scheme, besides possess- 
ing the desirable Delaunay properties, is that the procedure is very amenable 
to grid refinement and iterative node generation. The disadvantages of this 
particular scheme include the necessity to use high precision mathematics 
in certain situations and the need to introduce additional logic to preserve 
boundary integrity. The basic steps in the method are described below: 

1. Generate a fist of vertex locations. This list may include both bound- 
ary and interior vertices. In the present application of the scheme, 
an initial mesh is generated using only the boundary vertices and the 
final interior mesh is produced by using an iterative node generation 
scheme, similar to that suggested by Holmes and Snyder [HS88]. The 
iterative refinement scheme interrogates the existing mesh and locates 
points to improve the resolution, smoothness and aspect ratios of the 
cells. 

2. Create an initial temporary or virtual mesh that encloses the entire 
domain of interest. The present approach uses the triangulation of a 
quadrilateral (2D with 4 vertices) or hexahedron (3D with 8 vertices). 

3. Initially, introduce a vertex into the virtual mesh; and subsequently, 
into the existing triangulation from previous operations. 
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Adequately resolving the salient features of the solution field requires appro- 
priately clustering of the vertices. A higher density of grid points is required 
in regions where the geometry or solution vector is more complex. Although 
solution-adaptive grid refinement relaxes this requirement, it is still neces- 
sary to resolve the basic features of the computational field since the adaption 
parameters are computed from this initial field. Unlike techniques that use 
mappings to generate the mesh, unstructured grid generation techniques per- 
mit a different number of points on each of its boundaries. Therefore, the 
user can produce very dense vertex distributions on one boundary and ex- 
tremely coarse distributions on opposing boundaries. In addition, the regions 
may have numerous sides and even include embedded regions. 

Smooth variations in the sizes of the cells increases the accuracy of the nu- 
merical analysis. In particular, most numerical analysis techniques discretely 
represent the various terms in the mathematical expressions being solved by 
approximations whose accuracy degrades in the presence of rapid fluctuations 
in cell size. 

The aspect ratios and skewnesses of the cells can influence the accuracy 
and convergence characteristics of the numerical integration technique. The 
aspect ratios for triangles or tetrahedrons in the present study are defined 
by comparing the circumscribing radii with that of an optimal cell. An 
equilateral cell would have 0 skewness and a totally flat cell would have 
a skewness value of 1. A highly skewed cell decreases accuracy and slows 
convergence. 

In addition to the accuracy and convergence requirements of the numerical 
analysis technique, the robustness of the interior mesh generation technique is 
influenced by the quality of the boundary mesh. Smoothness is particularly 
critical in regions where surfaces intersect or where surfaces are in close 
proximity. For instance, the faces defining the bottom of a vehicle should 
not have dramatically different sizes than the faces defining the ground. 

Finally, since the purpose of the grid is to discretize the domain for numerical 
analysis, it is difficult to define the mesh quality requirements without some 
knowledge of the solution field. For example, highly skewed cells are tolerable 
in regions that contain no significant gradients. Nevertheless, the prime 
objective of the grid generation task is to cluster vertices in order to resolve 
salient features of the problem while producing smooth variations of low- 
aspect-ratio cells. 
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a decrease in node density on another area. The area from which nodes 
are moved usually has to be relatively near to the area into which they are 
moved (for not-trivial geometries), and this is not the place you tend to want 
to remove nodes from; you’d usually like to take them from somewhere far 
from the “action”. 

The obvious advantage of an unstructured mesh for solution adaption is 
that the node connectivity can change. This makes local node insertion 
and deletion possible. Less restricted node movement is also made possible 
because the node connectivity can be changed if a cell becomes too skewed. 

3.3.1 Point Insertion 

Rampant provides point insertion by allowing new points to be added at the 
midpoint of any edge. An edge is simply the line segment connecting two 
nodes. 

In two dimensions, an (interior) edge has a triangle on either side of it. 
When a new node is added at the mid-point of the edge, the original edge is 
removed and four new edges are added; two that connect the new node to 
the endpoints of the original edge, and two that connect the new node to the 
opposite vertices of the two adjacent triangular cells. Nodes are added on 
exterior boundaries in the same way, except that there is only one triangle to 
connect to (i.e., only three new edges). Periodic boundaries are logically the 
same as interior edges, but the implementation is more complicated because 
of the “cut” in computational space at the boundary; the insertion algorithm 
must wrap across the cut. 

In three dimensions, an edge has a “ring” of tetrahedra that surround it. As 
in 2D, when a new point is added along the central edge, each of the cells 
attached to it are divided into two, and the edges are reconnected. In 3D, 
the implementation becomes more complicated due to the variable number 
of cells attached to an edge and the intricate ways in which these rings of 
cells intersect solid boundaries, wrap across rotational periodic boundaries, 
and combinations thereof. 

Although any edge of a triangle or tetrahedron can be split in this manner, 
care must be taken in choosing the edge to be split so as not to create skewed 
cells. This is accomplished in Rampant by requiring the edge that is split 
to be the longest edge in each of the cells affected. If there is any longer 
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3.3.4 Solution-Adaptive Control 


Solution-adaptive mesh refinement is provided simply by controlling the 
point insertion and deletion procedures with information from the computed 
flow field. An adaption criterion is defined and applied to each of the cells 
in the domain. Depending on the value of the adaption function in the cell, 
the cell is marked for refinement, coarsening, or no change. A second pass 
is made through the domain, and in any cell marked for refinement, a new 
node is added at the mid-point of the longest edge in the cell. This may 
require refining neighboring cells first if the longest edge in the current cell 
is not also the longest edge in all the cells that contain it. After all cells 
marked for refinement have been refined, a third pass is taken through the 
domain and an attempt is made to coarsen any cell marked for coarsening. 
As described above, this is only implemented in 2D, and can only happen for 
certain cell configurations. In practice, the inability of the code to perform 
arbitrary coarsening is not much of a problem for steady solutions. It would 
be a significant problem for unsteady solutions where, for example, a shock 
might move across the domain leaving behind it a trail of small elements. 

There are two forms of adaption criterion available in Rampant. The first 
is one based on flow gradients. The adaption function is generated by first 
choosing a flow quantity such as pressure, density, Mach number, entropy, 
etc. Then node values of this function are computed by taking the average 
of the cell values in each of the cells that contain that node as a vertex. The 
value of the adaption criterion in each cell is then computed by subtracting 
the cell value of the function from the average of the three (or four in 3D) node 
values in that cell. This essentially results in an undivided Laplacian of the 
flow quantity. Two thresholds are then prescribed; any cell whose adaption 
value (magnitude) is below the lower threshold is marked for coarsening, 
and any cell whose adaption value is above the upper threshold is marked 
for refinement. 

The second adaption criteria is available for turbulent flows at walls and is 
based on the value of y + in the cell. Again, two thresholds are specified; 
any cell whose y + is less than the lower threshold is marked for coarsening, 
and any cell whose y + is greater than the upper threshold is marked for 
refinement. This provide a very easy way of generating a mesh with the 
proper spacing at walls. For the wall functions used in the solver, the first 
cell should have a y + less than 500 for the log-law to apply. In order to 
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Chapter 4 


Results 


4.1 Turbine Vane Cascade 

The first test case is a two-dimensional cascade of core turbine vanes. The 
three-dimensional annular geometry and experimental flow conditions for 
this case are described by Goldman and Seasholtz [GS82]. 

Figure 1 demonstrates the two-dimensional unstructured grid generation pro- 
cedure. First the geometry was defined in the PreBFC geometry modeler. 
The geometry consists of one curve for the inlet, one curve for the outlet, 
one curve for the periodic boundary, and three curves for the blade. Since 
finer spacing was desired at the trailing edge (and later on the upper surface), 
the blade was defined with three curves — upper surface, lower surface, and 
trailing edge — as a way to provide some simple grid spacing control. A single 
curve with non-uniform node spacing parameters could have been used (as 
is typical for a structured grid), but that degree of control was not opted for 
in this case. The disparities in the node spacing where the curves meet are 
easily smoothed away with a command in PreBFC that simply adds a few 
points to the end of the curve with the larger spacing. 

The three panels show the grid in its initial state, an intermediate state, 
and the final state. The initial grid is formed by tessellating the initial 
boundary points. This grid is a valid (constrained) Delaunay triangulation of 
the nodes; it simply has cells with very poor quality. To improve the quality 
of the elements, new points are inserted into the grid until the prescribed 
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4.2 Annular Turbine Cascade 


Next, a 3D inviscid computation was performed for the same geometry. The 
three-dimensional annular cascade of core turbine vanes is described by Gold- 
man and Seasholtz [GS82]. A very coarse mesh was used: 


9303 

cells 

19805 

faces 

2497 

nodes 

16763 

interior faces 

200 

outlet faces 

644 

periodic faces 

516 

wall faces (blade) 

671 

wall faces (hub) 

813 

wall faces (shroud) 

198 

inlet faces 


The memory required for the Euler solver was 184 bytes/cell, 28 bytes/face, 
and 24 bytes/node, for a combined cell-based requirement of 250 bytes per 
cell. 

Contours of Mach number for the inviscid solution are shown in Figures 12 
and 13. As expected, the results correspond directly to the previous 2D 
calculation. 


4.3 Laminar Flat Plate 

In order to ex ami ne the accuracy of the viscous stresses computed on the 
triangular grids, a laminar flat plate calculation was performed and compared 
to the Blasius similarity solution. Figure 14 shows the results for data taken 
at 5 locations along the plate. 

The plate starts at x=0 and extends to the outlet at x=10. The exact 
(Blasius) solution’s drag for the entire plate is 66.6 and the computed value 
was 64.4. The pressure residual converged 3 orders of magnitude in roughly 
3500 iterations. The grid had 16608 cells, 25425 faces, and 8818 nodes. 

The second panel shows all five data lines (the first two lines are indistin- 
guishable at that scale), and the fourth panel shows a close-up of just the first 
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52333 

cells 

109938 

faces 

12493 

nodes 

974 

inlet faces 

1862 

wall faces (shroud) 

1048 

wall faces (hub) 

5819 

wall faces (blade) 

2587 

periodic faces 

841 

outlet faces 

96807 

interior faces 

A (partially converged) solution was obtained on this starting mesh, and 
then the mesh was iteratively refined at solid boundaries as the solution was 

advanced until the value of in each cell was less than 500. The resulting 

mesh consisted of the following: 

180907 

cells 

370672 

faces 

37878 

nodes 

1112 

inlet faces 

4370 

wall faces (shroud) 

1134 

wall faces (hub) 

10249 

wall faces (blade) 

2894 

periodic faces 

851 

outlet faces 

350062 

interior faces 


Each iteration took 249 seconds on an SGI 4D/340VGX workstation. Ap- 
proximately 1000 iterations would be required to converge the solution on 
this mesh if the flow were reset to uniform conditions. The memory usage 
was 240 bytes/cell, 28 bytes /face, and 100 bytes/node for a combined cell- 
based requirement of 318 bytes per cell and a total of 58 Mbytes of grid 
and solution storage. The program code, graphics, interface, etc. add an 
overhead of about 6 Mbytes to the size of the process. 

The inlet stagnation properties were set at 
po = 101325 N/m2 
po = 1.226 kg/m3 

and the outlet static pressure at the hub was set at 
Phub = 101558 Pa = 1.0023 * p 0 
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of radial position. The open circles are experimental values, and the filled 
circles are computational values. Computed total pressure and temperature 
are consistently a few percent lower than experimental values, and the com- 
puted flow angle is consistently higher than the experimental values. This 
consistent behavior indicates that the discrepancies are due in large part to 
the different effective operating conditions for the two cases caused by the 
reduced blockage effects of the computed result. 
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